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ABSTRACT 

We present a comprehensive analysis of the ability of current stellar population 
models to reproduce the optical {ugriz) and near infra-red (JHK) colours of a small 
sample of well-studied nearby elliptical and SO galaxies. We find broad agreement be- 
tween the ages and metallicities derived using different population models, although 
different models show different systematic deviations from the measured broad-band 
fluxes. Although it is possible to constrain Simple Stellar Population models to a well 
defined area in age-metallicity space, there is a clear degeneracy between these param- 
eters even with such a full range of precise colours. The precision to which age and 
metallicity can be determined independently, using only broad band photometry with 
realistic errors, is A[Fe/H] ~ 0.18 and Alog Age ~ 0.25. To constrain the populations 
and therefore the star formation history further it will be necessary to combine broad- 
band optical-IR photometry with either spectral line indices, or else photometry at 
wavelengths outside of this range. 

Key words: galaxies: stellar content — galaxies: elliptical and lenticular, cD 



1 INTRODUCTION 

Knowledge of the star formation history of galaxies is key 
to our understanding of their formation and of environ- 
mental influences upon their properties. However, the re- 
lationship between observable quantities and the fundamen- 
tal parameters of the stellar population is complex. Observ- 
ables are usually interpreted by comparison with single-age, 
single-metallicity population (Simple Stellar Population or 
SSP) models. An SSP model is derived from a set of stel- 
lar tracks and isochrones, which are populated using an as- 
sumed Initial Mass Function (IMF). Empirical or theoretical 
spectra are summed along the isochrone, with appropriate 
weighting, producing the SSP spectral energy distribution. 
These Isochrone Synthesis Models were refined by Chariot 
& Bruzual (1991), and Bruzual & Chariot (1993) and found 
broadly to reproduce both spectroscopic and photometric 
properties of galaxies. However in detail the predictions are 
affected in a complex way by the population age, metallicity, 
abundance ratios such as the degree of enhancement of a- 
elements relative to solar ratios, and the IMF. Interpretation 
of gradients within individual galaxies, or trends with global 
parameters such as galaxy luminosity amongst samples, for 



instance in clusters, are subject to various degeneracies of 
which the best known is that between age and metallicity for 
old populations (Worthey 1994; Kodama & Arimoto 1997). 

Most attempts to break these degeneracies, and es- 
tablish luminosity-weighted stellar ages, metallicities and 
abundance ratios have involved obtaining high quality spec- 
tra of individual galaxies. A common technique is to com- 
pare plots of pairs of line indices with model grids (e.g. 
Trager et al. 2000a,b, 2008; Kuntschner et al. 2001; Pog- 
gianti et al. 2001a,b; Nelan et al. 2005; Sanchez-Blazquez et 
al. 2006). Other authors define combinations of indices de- 
signed to measure specific parameters (Vazdekis & Arimoto 
1999), or apply Principal Component Analysis or other non- 
parametric techniques to extract specific physical parame- 
ters from galaxy spectra (Panter et al. 2003, 2007; Tojeiro 
et al. 2007; Ocvirk et al. 2006a,b; Koleva et al. 2008). An 
alternative which represents a compromise between analy- 
sis of high signal-to-noise spectra, which can necessarily be 
applied only to limited samples, and broad band photome- 
try, is the technique of narrowband continuum photometry 
(Rakos et al. 2007, 2008). 

These techniques are limited to bright galaxies, where 
spectra or narrow band images of sufficiently high signal-to- 
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noise ratio can be obtained. For fainter samples, of dwarf 
or high redshift galaxies, or of star clusters associated with 
brighter galaxies, it is important to establish how precisely 
stellar population parameters can be measured from broad- 
band photometry alone. Optical broad-band photometry 
alone cannot break the age-metallicity degeneracy (e.g. Ko- 
dama and Arimoto 1999), however extending the spectral 
range appears to provide additional constraints. This has 
been demonstrated for optical - near infra-red colours by 
several studies, starting with Bothun et al. (1984) who ar- 
gued that J — K depends principally on metallicity, while 
B — H depends also on mean stellar age. The same colour 
combination was used used by Bothun & Gregg (1990) in 
an investigation of the star formation history of SO galaxies. 
Peletier, Valentijn and Jameson (1990) used U — B, B — V 
and V - K colours; Peletier & Balcells (1996) U - R vs. 
R- K; Bell & de Jong (2000) B - R vs R - K; and James 
et al. (2006) B — K vs. J — K in further studies of the 
metallicities and star formation histories of unresolved stel- 
lar populations in galaxies. Some success has also been found 
in the use of optical-ultraviolet colours, with Kaviraj et al. 
(2007a,b) combining SDSS optical and GALEX ultraviolet 
photometry in a study of the star formation histories of low- 
redshift early- type galaxies. 

One of the problems with all of these techniques is the 
model-dependence of the results. Each of the studies dis- 
cussed above will use a specific set of stellar tracks and 
isochrones, and a specific set of theoretical or empirical spec- 
tra. Differences between the results of these studies are of- 
ten the result of differences between models rather than 
between observations. In this paper we use a variety of 
popular SSP models to derive stellar population parame- 
ters from accurate photometry in (ugrizj H K) passbands 
of a sample of nearby E/SO galaxies. By using well cali- 
brated public datasets (SDSS and 2MASS) and photometry 
in consistent, well-defined apertures, the observational error 
is reduced and we can investigate in detail the uncertainty 
caused by the modelling process, the differences between 
popular sets of models, the difference between scaled-solar 
and Q-enhanced models, and that caused by different as- 
sumptions about the foreground extinction. Moreover for 
some galaxies our derived parameters can be compared di- 
rectly with spectroscopic studies (e.g. Trager et al. 2008). 

A somewhat similar study has been carried out by Em- 
inian et al. (2008). These authors use public catalogue pho- 
tometry in optical and near-infrared bands for a sample 
drawn from the SDSS spectroscopic catalog. Their study 
contains a high proportion of currently star forming galaxies, 
and their focus is on the correlation between star formation, 
dust, and near-IR colours, rather than the age and metal- 
licity of old populations. Their study relies on magnitudes 
from two different catalogue sources with different effective 
apertures, so while they use colours internal to each of the 
datasets (e.g. g — r and Y — K) they cannot use the SED 
over the entire range to constrain the population. The much 
lower precision of their photometry of individual galaxies is 
balanced by the much larger sample that they use. 

Conroy et al. (2009) use a new set of SSP models, cur- 
rently unavailable to us, to model the broad band SDSS 
and 2MASS catalogue magnitudes of samples of galaxies at 
low redshift and at redshift ~2. They emphasise the im- 
portance of understanding the Thermally Pulsing Asymp- 



totic Giant Branch (TP-AGB) phase of stellar evolution and 
of the IMF, and examine the effect of both on broad-band 
colours. They find that Horizontal Branch morphology im- 
pacts the derived population parameters less than either the 
TP-AGB and the IMF. They also conclude that the evolu- 
tion of a multi-metallicity population, particularly in pass- 
bands at wavelengths longer than V , is equivalent to that of 
a single-metallicity population whose metallicity is the mean 
of the multi-metallicity population. The study of Conroy et 
al. (2009), as with that of Eminian et al. (2008), relies on 
large samples of galaxies with individually much lower pho- 
tometric precision than that of our measurements. 

Lee et al. (2007), in a study which focuses on the com- 
posite populations found in spiral galaxies, find that because 
the youngest population (< 2 Gyr) dominates the light, 
broadband colours can partially break the age-metallicity 
degeneracy. However at older ages they find large scatter in 
the results obtained from models by different authors. 

A complementary approach is to apply similar tech- 
niques to extragalactic star clusters, in particular to the 
globular cluster populations of galaxies. This has the advan- 
tage that, although the cluster populations associated with 
galaxies can have multiple formation epochs and particu- 
larly metallicities, an individual cluster is more likely than a 
galaxy to contain stars of only a single age and metallicity. 
This is balanced against the disadvantage of using star clus- 
ters, that a comparatively small number of bright stars (on 
the AGB or on the upper main sequence, depending upon 
the age of the cluster) will contribute substantially to the lu- 
minosity, and that statistical fiuctuations in these numbers 
will cause fiuctuations in the photometry, and thus lead to 
additional uncertainty in the derived parameters of the stel- 
lar population. This is not a problem for the galaxies in our 
study, which will have typically 10* times as many stars in 
our defined apertures. 

Puzia et al. (2002), and Hempel & Kissler-Patig (2004a) 
applied V — I vs. V — K colour-colour plots to the study 
of unresolved stellar populations in globular clusters, in or- 
der to distinguish the metallicities of the separate cluster 
populations associated with the parent galaxies. Hempel & 
Kissler-Patig (2004b) add U band observations and show 
that the age resolution is enhanced, enabling separate forma- 
tion epochs to be distinguished as well. Kundu et al. (2005) 
find a substantial intermediate age population associated 
with the Virgo elliptical galaxy NGC4365, using glH pho- 
tometry alone. 

Anders et al. (2004) examine the ability of broad band 
(UBVRIJH) photometry to recover the age and metallic- 
ity of artificial clusters generated with GALEV SSP models 
(Schulz et al. 2002). They set observational errors to be 0.1 
mag in all passbands in their study, which is a somewhat 
larger value than we have in our observational data. They 
find that the U and B bands are critical for determining 
ages, and that the maximum available wavelength range is 
also important. Applying the same fitting technique to a 
sample of real clusters, de Grijs et al. (2005) find that for 
young star clusters, and now using HST filters and pass- 
bands, the estimates for the cluster age have an approxi- 
mately Gaussian distribution with <yiog{t) ^ 0.35. They do 
not attempt to constrain the metallicity of their clusters, 
but note that the effects of uncertainties in metallicity and 
extinction on the age determined are small. 
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Most of the Globular Cluster studies also rely on only 
one set of stellar population models, usually those of Bruzual 
& Chariot (2003, hereafter BC03). However Hempel et al. 
(2005), in their appendix, compare the model predictions for 
{V - I) and {V - H) for BC03 models and those of Maras- 
ton (2005), and find that the predictions for these colours 
do not differ significantly except for metal poor populations 
with age ^ 2Cyr. Pessev et al. (2008) present 2MASS-based 
JHK photometry of a sample of Magellanic Cloud globular 
clusters, and combine this with UBV photometry from the 
literature to undertake a comparison of the performance of 
four sets of stellar population models at distinguishing age 
and metallicity for this sample. They find that all models re- 
produce the colours for old clusters quite well, but that for 
young ages there are more substantial differences. Overall 
they conclude that the models of BC03 give the best quan- 
titative match to their data. This study is carried out at 
lower metallicity than ours, but comparison with our con- 
clusions will be interesting nevertheless. 



2 OBSERVATIONAL DATA 

Our sources of observational data are the image servers pro- 
vided by data release 6 of the Sloan Digital Sky Survey 
(SDSS; Adelman-McCarthy et al. 2008) and by the 2 Mi- 
cron All-Sky Survey (2MASS; Skrutskie et al. 2006). We 
have chosen a sample of bright galaxies in the Coma clus- 
ter, and a sample at intermediate luminosity from the Virgo 
cluster to test our techniques. These samples are limited 
by practical issues: for fainter galaxies in Coma the photon 
and read noise, which are the dominant noise source in our 
2MASS magnitudes, adversely affect the precision, while for 
brighter galaxies in Virgo the galaxy image is so large that 
it is impossible to determine sky levels on the size of image 
provided by the SDSS image server. To sample a range of 
absolute magnitude it was necessary to select galaxies from 
both clusters. 

SDSS and 2MASS images of a number of galaxies were 
downloaded from the image servers, and a sample of galax- 
ies was chosen. In the Coma cluster, as many galaxies as 
possible in common with the sample of Trager et al. (2008) 
were selected. In Virgo, galaxies were chosen where possible 
to be in the sample of the ACS Virgo cluster survey (Cote et 
al. 2004) and the catalogue of Michard (2005), note however 
that the majority of Michard's sample are too large for this 
study. Galaxies were rejected from the sample if they fell too 
close to an image boundary, to a neighbouring galaxy, or to 
a superimposed star to allow definition of a clean circular 
aperture centred on the galaxy. The final sample contained 
14 galaxies, eight from Coma and six from Virgo. In Table [T] 
we present the basic properties of this sample; galaxy types, 
redshifts and the Schlegel extinctions are taken from the 
NASA Extragalactic Database (NED), and the absolute V 
band magnitudes are taken from the Third Reference Cata- 
logue (RC3, de Vaucouleurs et al. 1991), assuming distance 
moduli of 35.0 for Coma, and 31.1 for Virgo. 

For each galaxy we measured aperture magnitudes in a 
series of apertures using the iraf task phot. Sky levels were 
determined in a number of regions of the image chosen to be 
free of foreground stars. Galaxy centres were determined in 
each image using the centroiding function within phot. For 
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Table 1. Basic propert: 


ies of the sample galaxies. 




Galaxy 


Type 


Redshift 


Av 


My 




(NED) 




(Schlegel) 


(RC3) 


IC3501 


d:El 


0.0055 


0.092 


-17.37 


NGC4318 


E 


0.0041 


0.083 


-17.88 


NGC4515 


SO- 


0.0032 


0.103 


-18.60 


NGC4551 


E: 


0.0039 


0.128 


-19.21 


NGC4564 


E6 


0.0038 


0.116 


-20.10 


NGC4Sfi7 


E3 


0.0162 


0.034 


-20.57 


NGC4872 


SBO 


0.0241 


0.030 


-20.62 


NGC4871 


SABO/a 


0.0227 


0.030 


on Qn 
— zU.oy 


NGC4873 


SAO 


0.0194 


0.027 


-20.93 


NGC4473 


E5 


0.0075 


0.094 


-20.99 


NGC4881 


EO 


0.0225 


0.036 


-21.47 


NGC4839 


cD 


0.0246 


0.032 


-22.97 


NGC4874 


cD 


0.0241 


0.028 


-23.35 


NGC4889 


cD/E4 


0.0217 


0.032 


-23.54 



Notes to Table]!]; Column 1 gives the galaxy name, column 2 the 
type from NED. Column 3 gives the redshift, again from NED, 
column 4 the extinction Ay using the Schlegel (1998) maps. 
Column 5 is a y band absolute total magnitude for the galaxy 
from RC3, using the distance modulus estimates given in the 
text. The IC3501 My is determined from interpolation of g and 
r flux densities. 

each galaxy we then chose an optimum aperture, in the range 
20 - 40 arcsec diameter, within which to measure the mag- 
nitude in all passbands. All apertures are large enough that 
the different point spread functions of SDSS and 2MASS 
data do not have a significant effect. The optimum aperture 
chosen was in each case a compromise between the size of 
the galaxy, the proximity of neighbours, and the additional 
noise in the 2MASS magnitudes in larger apertures. 

Magnitude zeropoints were taken directly from the fits 
headers of the 2MASS images, however for the SDSS data 
the zero point keyword (FLUX20) in the fits header of the 
downloaded image files (fpC files) has not been calculated 
using the final calibration informatiorQ. Instead we recali- 
brated the fpC images using the zeropoint, extinction cor- 
rection, and airmass given in the tsField files for the relevant 
Run and Field number. This recalibration gives a substan- 
tial correction to the magnitudes calculated using the image 
headers, particularly in the u, g and z bands where the cor- 
rections are in the range 0.05 - 0.1 magnitudes. 

2.1 Homogenising galaxy photometry 

In order to compare the models to our galaxy observations, 
the measured magnitudes must be homogenised to the same 
photometric system, for which we choose the AB magnitude 
system. Our optical photometry is carried out in the SDSS 
DR6 system, however there are sky position dependent dif- 
ferences between DR6 and DR7 zero points (see Padman- 
abhan et al. 2008). We transform our photometry to the 
DR7 photometric system by adding the difference between 
zero points calculated for each galaxy from the difference be- 
tween DR6 and DR7 catalogue magnitudes. The SDSS DR7 
photometry is brought onto the AB magnitude system by 
subtracting 0.04 from the u and adding 0.02 to the z band 

^ http:/ /www. sdss.org/DR6/algorithms/fluxcal. html 
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Table 2. Conversions between 2MASS (Vega) and AB magni- 
tude systems for the JHK bands, taken from Cohen, Wheaton 
&; Megeath (2003). To apply the corrections, niAB = ^Vega + c, 
where c is the relevant correction from the table, shown to 3 deci- 
mal places. The uncertainties for each of the corrections are added 
in quadrature to the photometric errors in the J, H &i Kg bands. 



Band J H Ks 

Correction 0.894 1.374 1.840 
Uncertainty 0.022 0.022 0.022 



Magnitudes are not K-corrected, as the fitting proce- 
dure we use integrates the redshifted model flux over the 
rest wavelength filter passbands, however they are corrected 
for Galactic extinction, using the Ay estimates of Schlegel 
et al. (1998), as tabulated in NED, and converted to SDSS 
passbands using the prescription of Kim and Lee (2007), 
and to JHK using the Schlegel et al. prescription which 
is in turn derived from the functional form of Clayton & 
Cardelli (1988). The final magnitudes and errors used for 
the fitting are presented in Table 3. 



magnitudes as recommended on the SDSS web page fl We 
note that Eisenstein et al. (2006) , in their work on hot white 
dwarfs, recommend a slightly difi'erent AB transformation 
vector, the main difference of which is to add approximately 
0.015 to the i band magnitudes. This vector might improve 
marginally many of our fits, however as the SDSS project 
appears not to have adopted the Eisenstein et al. transfor- 
mation vector we do not feel justified in doing so. 

The Infrared 2MASS magnitudes make use of the Vega 
magnitude system; these are brought on to the AB system by 
adopting the conversions from Cohen, Wheaton & Megeath 
(2003), shown in table [2] As a check we also convolved the a 
Lyiee (Vega) spectral energy distribution of Kurucz (1993) 
with the 2MASS JHKs response curves and converted the 
results to AB magnitudes, in effect calculating the magni- 
tudes of Vega in the AB system. After taking into account 
the AB magnitude of Vega (0.03), our values are in agree- 
ment with those of Cohen, Wheaton & Megeath (2003). 



2.2 Photometric errors and corrections 

For each magnitude we calculated an error, allowing for 
three sources of uncertainty. The first was the error in the 
measurement process, which in turn results from errors in 
the determination of the sky level. The measurement error 
was estimated from the SDSS and 2MASS images by re- 
peated measurements of the magnitudes from the same im- 
age. In the 2MASS images this is a significant, albeit never 
dominant, source of error. The second source of error comes 
from noise in the image. This is calculated from the image 
gains given in the fits headers, and the image statistics as 
measured in the apertures in the case of SDSS passbands, 
and as given in the image headers in the case of 2MASS. For 
the 2MASS data the errors are corrected to allow for the re- 
sampling and smoothing applied to the data before delivery, 
using the prescription given on the 2MASS web pages|f] 

The third source of error is the error on the photo- 
metric zero point, and this is taken from the web pages of 
SDS^I and 2MAS^| respectively. For the 2MASS data we 
consider also the uncertainty in the conversion of 2MASS to 
AB magnitudes, as given in Table (2] These error sources are 
all added in quadrature. The zero point error always domi- 
nates for all SDSS passbands, whereas the image noise error 
dominates in most cases for 2MASS bands. 



^ http:/ /www.sdss.org/dr6/algorithms/fluxcal.html#sdss2ab 

^ http: / / www.ipac.caltech.edu/2mass / releases / allsky / doc / sec6_8i 

* http:/ /www. sdss.org/dr5/algorithms/fiuxcal. html 

^ http: / / www.ipac.caltech.edu/2mass /releases / allsky / doc / sec4_8 



3 STELLAR POPULATION MODELS 

We have chosen to test seven popular sets of stellar popu- 
lation models against the data presented in Section [2l All 
models produce a series of synthetic spectra at points on a 
grid in age-metallicity space which are compared with the 
broad-band photometry as described in Section |4] For con- 
venience, the basic parameters of each model are presented 
in Table 4, and we describe here briefiy the further relevant 
details, in particular the treatment of the TP-AGB stars 
which are generally not included in the adopted stellar evo- 
lutionary tracks. We refer the reader to the original papers 
describing these models for a full description of each model. 

(i) Bruzual & Chariot (2003) models (BC03) use the 
Padova 1994 evolutionary tracks (Alongi et al. 1993; Bressan 
et al. 1993; Fagotto et al. 1994a,b; Girardi et al. 1996) with 
the TP-AGB and Post Asymptotic Giant Branch (PAGB) 
phases treated according to Vassiliadis & Wood (1993, 
1994). We use the STELIB spectral library in our modelling, 
although our results do not change if we use the BaSeL-2.0 
library. Neither stellar library contains carbon stars or up- 
per TP-AGB stars and the spectra for these stars are con- 
structed from observational libraries as described by BC03. 

(ii) PEGASE models (Fioc & Rocca-Volmerange 1997, 
1999) also use stellar tracks from the Padova group, and 
pseudo-tracks for the TP-AGB phase as proposed by Groe- 
newegen & de Jong (1993). 

(iii) Starburst99 (Vazquez & Leitherer 2005) is of par- 
ticular interest in the study of young stellar populations, 
although it is designed to model populations of any age. It 
uses the Padova 1994 stellar tracks and incorporates TP- 
AGB stars in a similar manner to BC03, as described by 
Vazquez & Leitherer (2005). Spectra were computed using 
the PAULDRACH/HILLIER option, which uses the BaSeL-2.0 
atmospheres for all except the hottest stars. Starburst99 
models using Geneva high mass-loss isochrones were also 
computed, but provided a very poor fit. 

(iv) GALEV (Anders & Fritze-van Alvensleben 2003) is 
an evolutionary synthesis code written to study the spectral 
and chemical evolution of galaxies, underpinned by a set of 
SSP models (Schulz et al. 2002; Anders et al. 2004) using 
Padova isochrones and BaSeL-2.0 atmospheres. 

(v) The SPEED models are described by Jimenez et al. 
(1995, 2004) and adopt a stellar evolution code based upon 
Eggleton (1971, 1972) but with improvements, particularly 
in the modelling of the HB and AGB. Mass loss at the RGB 
tip is modelled using a distribution of values of the Reimers 

ihtm{19^5) mass loss parameter r], and this distribution is cal- 
ibrated against the observed HB morphologies of Galactic 
html star clusters (Jimenez et al. 1996). The AGB is modelled 
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Table 3. Photometry of the sample galaxies. 



Galaxy 


D 

arcsec 


n° 






a° 








.0 






i° 








t° 




(Tz 


J° 




H° 




K° 




IC3501 


20 


16. 


,460 


0, 


,033 


14, 


.991 


0, 


.020 


14, 


,318 


0, 


,020 


13, 


.995 





.020 


13, 


,795 


0, 


.030 


13.493 


0, 


,027 


13.332 





.047 


13.577 


0, 


.027 


NGC4318 


30 


15. 


,393 


0, 


,031 


13, 


.824 





.020 


13, 


122 


0, 


,020 


12, 


.756 





.020 


12, 


,557 





.030 


12.248 


0, 


,025 


12.045 





.027 


12.313 


0, 


.030 


NGC4515 


30 


14. 


,914 


0, 


.030 


13, 


.394 


0, 


.020 


12, 


,714 


0, 


,020 


12, 


.355 





.020 


12, 


,163 





.030 


11.854 


0, 


.023 


11.704 


0, 


.026 


11.976 


0, 


.024 


NGC4551 


40 


14. 


,521 


0, 


.031 


12, 


.798 


0, 


.020 


12, 


,049 


0, 


,020 


11, 


.660 





.020 


11, 


,416 





.030 


11.008 


0, 


,022 


10.822 


0, 


.023 


11.063 


0, 


.022 


NGC4564 


40 


13. 


,845 


0, 


.030 


12, 


.111 


0, 


.020 


11, 


,345 


0, 


,020 


10, 


.942 





.020 


10, 


,694 





.030 


10.278 


0, 


,022 


10.052 





.023 


10.292 


0, 


.022 


NGC4867 


20 


16. 


,855 


0, 


.034 


15, 


.150 





.021 


14, 


,363 


0, 


,020 


13, 


.967 





.020 


13, 


,696 





.031 


13.297 


0, 


,027 


13.035 





.030 


13.256 


0, 


.031 


NGC4872 


20 


16. 


,780 


0, 


.033 


15, 


.070 





.020 


14, 


,267 


0, 


,020 


13, 


.865 





.020 


13, 


,596 





.030 


13.210 


0, 


,026 


12.963 





.029 


13.164 


0, 


.030 


NGC4871 


20 


16. 


,868 


0, 


.033 


15, 


.105 


0, 


.020 


14, 


,294 


0, 


,020 


13, 


.893 





.020 


13, 


,620 


0, 


.030 


13.233 


0, 


,026 


12.964 


0, 


.029 


13.190 


0, 


.030 


NGC4873 


20 


16. 


,802 


0, 


.033 


15, 


.135 


0, 


.020 


14, 


,354 


0, 


,020 


13, 


.955 





.020 


13, 


,685 


0, 


.030 


13.351 


0, 


,027 


13.120 





.031 


13.332 


0, 


.032 


NGC4473 


40 


13. 


,268 


0, 


.030 


11, 


.509 





.020 


10, 


,737 


0, 


,020 


10, 


.335 





.020 


10, 


,083 





.030 


9.666 


0, 


.022 


9.466 





.022 


9.676 


0, 


.021 


NGC4881 


30 


16. 


,255 


0, 


.034 


14, 


.527 


0, 


.020 


13, 


,716 


0, 


,020 


13, 


.314 





.020 


13, 


,039 


0, 


.030 


12.634 


0, 


.025 


12.399 





.027 


12.605 


0, 


.028 


NGC4839 


40 


15. 


,501 


0, 


.032 


13, 


.695 


0, 


.020 


12, 


,877 


0, 


,020 


12, 


.460 





.020 


12, 


,177 


0, 


.030 


11.739 


0, 


,023 


11.502 


0, 


.024 


11.703 


0, 


.024 


NGC4874 


40 


15. 


,353 


0, 


.031 


13, 


.565 


0, 


.020 


12, 


,741 


0, 


,020 


12, 


.317 





.020 


12, 


,045 


0, 


.030 


11.633 


0, 


.023 


11.390 


0, 


.024 


11.582 


0, 


.024 


NGC4889 


40 


14. 


,889 


0, 


.031 


13, 


.054 


0, 


.020 


12, 


,220 


0, 


.020 


11, 


.789 





.020 


11, 


,509 


0, 


.030 


11.061 


0, 


.023 


10.819 





.022 


10.997 


0, 


.022 



Notes to Table S: Column 1 gives the galaxy name, column 2 the aperture diameter for the photometry in arcseconds. Columns 3 to 18 give the extinction 
corrected magnitude within that aperture, and the formal error on that magnitude, in ugrizJH K passbands in that order. All magnitudes are on the AB system. 



Table 4. Key parameters of the stellar population synthesis models used. Additional details, including the different treatments of post-main sequence evolution, are given in the text. 



Name 



Reference 



Spectral 
library ■'■ 



Range of Z 
(number) 



Range of ages 
(Gyr) (number) 



Availability of models 



BC03 Bruzual & Chariot (2003) Chabrier (2003) STELIB 0.0001-0.05 (6) 0.0001-20 (221) 

PEGASE-2.0 Fioc & Rocca-Volmerange (1997,1999) Kroupa (2001) BaSeL-2.0 0.0001-0.1 (7) 0.5-13 (25) 

Starburst99 Vazquez & Leitherer (2005) Kroupa (2001) BaSeL-2.0 0.0004-0.05 (5) 0.00001-15 (100) 

GALEV Anders & Fritze-van Alvensleben (2003) Kroupa (2001) BaSeL-2.0 0.0004-0.05 (5) 0.004-16 (51) 

SPEED Jimenez et al. (1995,2004) Kroupa (2001) Kurucz (1993) 0.001-0.05 (4) 0.001-14 (40) 

BaSTI Pietrinferni et al. (2004) Kroupa (2001) see text 0.0003-0.04 (10) 0.03-15 (55) 

Maraston Maraston (2005) Kroupa (2001) BaSeL-2.0 0.0001-0.07 (6) 0.001-15 (40) 



http : //www. cida. ve/~bruzual/bc2003 
http : //www2 . lap . f r/users/f ioc/PEGASE. html 
http : //www. stsci . edu/science/starburst99 
http : //www . galev . org 

http: //www. ice . csic . es/personal/ j imenez/SPEED/ 
http: //albione . oa-teramo . inaf . it/ 
http : //www. dsg. port . ac .uk/~marastonc/ 



^The STELIB hbrary is described by Le Borgne et al. (2003) and the BaSeL-2.0 library by Lejeune et al. (1997, 1998) 
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as described by Jorgensen (1991). Spectra are largely taken 
from Kurucz (1993), but the authors compute their own LTE 
stellar atmosphere models for stars cooler than 4000 K. 

(vi) Maraston (2005) presents a set of models determined 
from an evolutionary synthesis code (Maraston 1998) based 
upon the tracks of Cassisi et al. (1997a, b; 2000) for main 
sequence stars, and upon the Fuel Consumption Theorem 
(Renzini & Buzzoni 1986) for post-MS evolution, calibrated 
against clusters in the Magellanic Clouds. 

The BaSeL2.0 library of stellar spectra is used, with TP- 
AGB stars taken from Langon & Mouhcine (2002). Maras- 
ton's models are available with two horizontal branch mor- 
phologies (RHB and BHB) and we primarily consider the 
RHB models which are more appropriate for metal-rich pop- 
ulations. However, the existence of a blue HB morphology 
in a small number of metal-rich Galactic globular clusters 
(Busso et al. 2007) leads us to consider the BHB models as 
well. 

Maraston et al. (2008) discuss a set of related models 
which use the Pickles (1998) stellar library instead of BaSeL- 
2.0, but these were not available to us at the time of writing. 

(vii) The BaSTI (Bag of Stellar Tracks and Isochrones) 
models are based on the work of Pietrinferni et al. (2004), 
extended to cover TP-AGB stars by Cordier et al. (2007). 
Model spectra are constructed from the BaSTI isochrones 
as described by Percival et al. (2008), using the 'low reso- 
lution' version of the model spectra. These are based upon 
the model atmospheres of Castelli & Kurucz (2003), supple- 
mented by models from the BaSeL-3.2 library (Westera et 
al. 2002) for cool stars (T < 3500 K), and by the empirical 
spectra of Langon & Mouhcine (2002) for AGB carbon stars. 

Because these models are 'in-house', we are able to test 
different prescriptions for a-enhancement and HB morphol- 
ogy. In addition to the scaled-solar models (Pietrinferni et 
al. 2004) we analyse also the [a/Fe] — -1-0.4 models and 
isochrones of Pietrinferni et al. (2006) . We also analyse two 
different values (-1-0.2 and -1-0.4) of the Reimers (1975) mass- 
loss parameter, 77. Discussion of the effects of varying [a/Fe] 
and ri is presented in Section [5.21 



4 INTEGRATION AND FITTING 
PROCEDURES 

4.1 Calculating model fluxes 

In order to compare our observations with the results of 
SSP modelling, we convert from our observed AB magni- 
tudes and magnitude errors to fluxes. To extract compara- 
ble fluxes, Sn, from the SSP SEDs of each species there are 
a number of steps. First, we multiply the wavelength values 
of the given SSP model by a factor of {1 + z) to account 
for the redshift of the galaxy in question. Second we rebin 
linearly the transmission function of each filter, TA.n, onto 
the wavelength scale of the models. 

For each SSP SED /a, filter n, and filter transmission 
function Tx,n, we then integrate the total flux in the pass- 
band using trapezoidal summation: 



fxXTx-ndX 



XT\^„dX 



(1) 



Next, we convert the observed AB magnitudes to fx 
using: 



(2) 



where /o is the AB magnitude zero point of 3631 Jy, and 
Xeff,n is the effective, or pivot, wavelength of the passband, 
defined by: 



_jTxXdX 



(3) 



We are now in a position to compare our observations with 
the model SSP fluxes. 



4.2 X 



^ fitting 



In order to quantify the extent to which the galaxies in our 
sample are described by the array of models mentioned in 
section 3, we employed a simple reduced test. For a given 
galaxy with observed fluxes fi with each flux having an as- 
sociated error ai in n bands, the extent to which they are 
well described by a model with fluxes Si is approximated by 
the reduced x'^i calculated as follows; 



2 



in- 



^ (NSi - f^f 

i=0 



(4) 



where A*' represents the normalisation of the model magni- 
tudes and K represents the number of constraints associated 
with the models. The normalisation factor N is calculated 
by analytically solving the equation Ax^/AN = 0. The value 
of K used for these calculations is 3, due to the treatment of 
a galaxy model's age, metallicity, and flux density normali- 
sation as free parameters. Model SSPs that approximate our 
galaxy observations well will then have values of reduced x^ 
approaching unity, while for those models which do not ap- 
proximate our observations well the values will be higher, 
depending on the 'goodness of fit'. 

In order to discriminate between the large numbers of 
models for a given galaxy, the best fit value of x^ is calcu- 
lated for every SSP in each model species (e.g. BC03, M05, 
etc.). In this way we build up a database of best-fit x^ for 
all values of SSP age and metallicity for each galaxy in our 
sample. Making use of this database, we can then compare 
the results of our x^ fitting to those of other studies (e.g. 
Trager et al., 2008). The properties of the best fit models to 
each of the galaxies in our sample are shown in tables 5 and 
6. 



5 ANALYSIS 

Remembering always that we are observing potentially com- 
posite populations, and comparing them with SSP models, 
we examine in turn the model sets, and how well they fit to 
the data. 

(i) BC03 models. The BC03 models often provide the low- 
est values of a^nd thus objectively the best fit for the low 
luminosity Virgo galaxies in particular. In all cases this best 
fit comes from the [Fe/H] = 0.093 model, at an age between 
2.5 and 12 Gyr, with the massive Coma galaxies in general 
showing the oldest ages. The residuals around the fit show 
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Table 5. Best fit BC03, PEGASE, Starburst99, Galev, SPEED and Maraston models 



Model BC03 PEGASE Starburst99 GALEV SPEED Maraston RHB Maraston BHB 

xl Age [^] xl Age xl Age xl Age [^] xl Age [^] xl Age 



IC3501 


0, 


,55 


2.50 


0, 


,093 


2, 


.20 


3, 


,0 


-0.204 


1, 


,46 


2. 


,72 


-0.40 


0, 


,93 


3.16 


-0.40 


3.21 


5.5 


-0.70 


1, 


,40 


1, 


.5 


0, 


,00 


4, 


,06 


10, 


,0 


0, 


,00 


NGC4318 


1, 


,40 


3.00 


0, 


,093 


2, 


.33 


5, 


,0 


-0.204 


2, 


,86 


2. 


,72 


0.00 


1, 


,80 


6.31 


-0.40 


3.71 


9.0 


-0.70 


4, 


,65 


1, 


.5 


0, 


,67 


6, 


,81 


10, 


,0 


0, 


,00 


NGC4515 


2, 


,06 


2.75 


0, 


,093 


1, 


.93 


4, 


,0 


-0.204 


2. 


,60 


3. 


,37 


-0.40 


2, 


,72 


5.01 


-0.40 


6.96 


5.5 


-0.70 


2. 


,89 


2, 


.0 


0, 


,00 


3, 


,21 


10, 


,0 


0, 


,00 


NGC4551 


1, 


,56 


7.50 


0, 


,093 


2, 


.67 


3, 


,0 


0.176 


2. 


,40 


2. 


,92 


0.40 


4, 


,00 


6.31 


0.00 


12.67 


14.0 


-0.70 


3, 


,44 


2, 


.0 


0, 


,67 


4, 


,26 


10, 


,0 


0, 


,35 


NGC4564 


1, 


,89 


8.50 


0, 


,093 


2, 


.88 


3, 


,5 


0.176 


2, 


,84 


3. 


,14 


0.40 


2, 


,62 


6.31 


0.00 


23.00 


14.0 


-0.70 


3, 


,59 


2, 


.0 


0, 


,67 


5, 


,92 


15, 


,0 


0, 


,35 


NGC4867 


2, 


12 


8.00 


0, 


,093 


2, 


.67 


3, 


,5 


0.176 


2, 


,40 


6. 


,86 


0.00 


1, 


,79 


3.98 


0.00 


12.40 


14.0 


-0.70 


3, 


,30 


2, 


.0 


0, 


,67 


7, 


,21 


15, 


,0 


0, 


,35 


NGG4872 


1, 


,23 


8.00 


0, 


,093 


2, 


.25 


6, 


,0 


0.000 


2, 


,22 


6. 


,39 


0.00 


1, 


,99 


6.31 


0.00 


10.54 


14.0 


-0.70 


3, 


10 


2, 


.0 


0, 


,67 


7, 


,50 


10, 


,0 


0, 


,35 


NGG4871 


0, 


,94 


8.50 


0, 


,093 


2, 


.41 


8, 


,0 


0.000 


2, 


18 


6. 


,86 


0.00 


2, 


,74 


6.31 


0.00 


17.20 


14.0 


-0.70 


4, 


12 


2, 


.0 


0, 


,67 


5, 


,95 


15, 


,0 


0, 


,35 


NGG4873 


0, 


,83 


6.50 


0, 


,093 


1, 


.86 


4, 


,5 


0.000 


1, 


,87 


3. 


,89 


0.00 


4, 


,59 


13.18 


-0.40 


4.24 


14.0 


-0.70 


5, 


,62 


2, 


.0 


0, 


,67 


4, 


,99 


10, 


,0 


0, 


,35 


NGG4473 


1, 


,53 


8.50 


0, 


,093 


3, 


.40 


3, 


,5 


0.176 


2, 


,57 


3. 


,14 


0.40 


3, 


,30 


6.31 


0.00 


25.11 


14.0 


-0.70 


3, 


,42 


2, 


.0 


0, 


,67 


5, 


,48 


15, 


,0 


0, 


,35 


NGG4881 


1, 


13 


8.50 


0, 


,093 


2, 


.22 


3, 


,5 


0.176 


1, 


,90 


6. 


,86 


0.00 


2, 


,04 


6.31 


0.00 


15.59 


14.0 


-0.70 


3, 


,55 


2, 


.0 


0, 


,67 


6, 


,34 


15, 


,0 


0, 


,35 


NGG4839 


1, 


,76 


10.00 


0, 


,093 


2, 


.25 


4, 


,0 


0.176 


2, 


,40 


8. 


,49 


0.00 


2, 


,49 


7.94 


0.00 


31.70 


4.5 


0.00 


7, 


,83 


5, 


.0 


0, 


,35 


4, 


,81 


15, 


,0 


0, 


,35 


NGG4874 


1, 


,20 


9.50 


0, 


,093 


2, 


.54 


4, 


,5 


0.176 


2, 


,96 


8. 


,49 


0.00 


2, 


,77 


7.94 


0.00 


33.95 


2.5 


0.00 


7, 


,53 


5, 


.0 


0, 


,35 


5, 


.49 


15, 


,0 


0, 


,35 


NGG4889 


2, 


,25 


12.00 


0, 


,093 


3, 


.24 


6, 


,0 


0.176 


2, 


,99 


9. 


,12 


0.00 


3, 


,74 


10.00 


0.00 


29.21 


8.0 


0.00 


9, 


,68 


6, 


.0 


0, 


,35 


9, 


,58 


15, 


,0 


0, 


,35 



Notes to Table 5: For each model set, as described in text, there are three columns in which we tabulate the best fit xt, and the age (in Gyr) and [Fe/H] for the best fit model. 
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Table 3. BaSTI models for different a-enhancement and mass loss parameters 















BaSTI 














[a/Fe] 


= 0.0, T) 


= 0.2 


[a/Fe] 


= 0.0, r) 


= 0.4 




= 0.4, ri 


= 0.2 


[a/Fe] 


= 0.4, ri 


= 0.4 




xl 


Age 


[^1 


xl 


Age 


I H 1 




Age 


[^1 


xl 


Age 


[^] 


IC3051 


1.61 


4.50 


-0.25 


1.47 


4.50 


-0.25 


1.04 


8.50 


-0.60 


0.95 


14.00 


-0.60 


NGC4318 


1.65 


8.00 


-0.25 


1.58 


8.00 


-0.25 


1.56 


13.50 


-0.60 


2.27 


5.50 


-0.29 


NGC4515 


2.66 


6.00 


-0.25 


2.56 


6.00 


-0.25 


2.67 


10.00 


-0.60 


2.87 


10.50 


-0.60 


NGC4551 


2.62 


2.75 


0.40 


2.50 


2.75 


0.40 


2.76 


7.50 


-0.09 


2.69 


8.00 


-0.09 


NGC4564 


2.78 


3.00 


0.40 


2.71 


3.00 


0.40 


1.31 


6.50 


0.05 


1.40 


10.00 


-0.09 


NGC4S67 


2.91 


3.00 


0.40 


2.82 


3.00 


0.40 


0.54 


6.00 


0.05 


0.57 


6.00 


0.05 


NGC4872 


2.43 


3.00 


0.40 


2.37 


3.00 


0.40 


0.73 


8.50 


-0.09 


0.74 


9.00 


-0.09 


NGC4871 


2.09 


3.00 


0.40 


2.07 


3.00 


0.40 


0.89 


10.00 


-0.09 


0.86 


10.00 


-0.09 


NGC4873 


2.05 


3.00 


0.26 


2.05 


3.00 


0.26 


0.91 


11.00 


-0.29 


0.99 


11.00 


-0.29 


NGC4473 


2.37 


3.00 


0.40 


2.33 


3.00 


0.40 


2.36 


6.50 


0.05 


2.39 


10.00 


-0.09 


NGC4881 


2.60 


3.00 


0.40 


2.54 


3.00 


0.40 


0.77 


6.50 


0.05 


0.78 


10.00 


-0.09 


NGC4839 


3.90 


6.00 


0.40 


2.70 


5.00 


0.40 


1.17 


8.50 


0.05 


1.18 


8.50 


0.05 


NGC4874 


3.38 


4.50 


0.40 


2.20 


5.00 


0.40 


1.28 


11.50 


-0.09 


1.31 


12.00 


-0.09 


NGC4889 


4.46 


6.50 


0.40 


4.37 


6.50 


0.40 


1.34 


10.00 


0.05 


1.35 


10.50 


0.05 



Notes to Table \^ Columns 2 — 13 - Best fitting BaSTI models for two values of [a/Fe], and for two values of the Reimers (1975) mass 
loss parameter rj. For each combination of parameters, there are three columns in which wo tabulate the best fit xl^ ^^nd the age (in 

Gyr) and [Fe/H] for the best fit model. 



a systematic pattern, in particular (r — i) is predicted too 
blue by 0.01 - 0.03 magnitudes for the lower luminosity Virgo 
galaxies, and 0.03 - 0.05 magnitudes for the more massive 
Coma galaxies. BC03 models provide too coarse a metal- 
licity grid to investigate age-metallicity degeneracy, and in 
particular their one supersolar metallicity is [Fe/H] — 0.559, 
which is apparently too metal rich to fit any of our sample. 

In Figure [T] we show the dependence of xl on ^S^) the 
best-fit model spectrum, and the residuals between data and 
best fit model, all for NGC4881, one of the best fitting galax- 
ies. Similar plots for all galaxies and all model sets are avail- 
able in the online data accompanying this paper. 

(ii) PEGASE models. These models provide higher values 
of xl than BC03, usually because they underpredict the flux 
in the i band. PEGASE models predict {r — i) to be too blue, 
and {i — z) to be too red, both by as much as 0.1 magnitude. 
The solar metallicity models show a minimum at generally 
slightly younger ages than BC03, but there is quite often 
a better fit from an interpolated model at [Fe/H] = -0.204 
or -1-0.176. The xl curves for these models show clearly the 
age-metallicity degeneracy, where the supersolar model gives 
a fit a factor 1.5 - 2 younger, and the subsolar one a factor 
of ^2 older. 

(iii) Starburst99 models. Using Padova isochrones, these 
models give very similar results to PEGASE, in that they 
predict the same (r — i) and (i — z) discrepancy, they show 
very similar ages, and the same age-metallicity degeneracy. 
They have less metallicity resolution than PEGASE, and 
sometimes this results in a higher best value of xl- Star- 
burst99 models using Geneva isochrones fit substantially less 
well in all cases, and we do not consider these models fur- 
ther. 

(iv) GALEV models. GALEV in general does not give 
quite as good a fit as PEGASE, but the best fits usually 
occur at a similar age and metallicity. GALEV ages are sim- 
ilar to PEGASE ages, and in general older than Starburst99 
ages. These models again predict both (r — i) to be too blue 
and (i — z) to be too red, by 0.05 - 0.1 magnitude. 

(v) SPEED models. The best fit SPEED models are al- 



most always old (often at their limit of 14 Gyr), and metal- 
poor (Z=0.001 or 0.004). They also overpredict the u band 
flux by 0.1 mag or more, and as with other model sets pre- 
dict (r - i) to be too blue by 0.05 - 0.1 mag. The xl values 
are consequently high compared with other models. 

(vi) Maraston (2005) models. The RHB models provide a 
less good fit to the broad band data than BC03, PEGASE, 
Starburst99 or BaSTI, with < 3 for only two of the four- 
teen galaxies. The general pattern of the residuals is that 
(r — i) is predicted too blue, by 0.15 magnitude. For the 
Coma galaxies (J — K) is also too blue, by around 0.05 
magnitude. The best fit ages are also uncomfortably young, 
in the range 1.5 - 3 Gyr even for massive cluster ellipticals. 
Such ages are in conflict with those derived from line indices 
(Trager et al. 2008). The BHB models are only appropriate 
for ages > 10 Gyr and would generally be expected to fit low 
metallicity populations, so it is somewhat surprising that the 
[Fe/H] — 0.35 BHB model provides a better fit to the data 
than any of the RHB models, for the massive Coma galaxies 
in particular. 

(vii) BaSTI models. The standard scaled-solar abundance 
models give slightly higher values of xl than BC03, and 
similar to PEGASE or Starburst99. The residuals show an 
overprediction of the r band fiux, of 0.03 - 0.08 magnitudes. 
As with PEGASE, BaSTI gives younger ages than BC03, 
largely because it has two supersolar metallicities ([Fe/H] — 
+0.26 and -1-0.40). In section 15.21 we discuss the effect of 
making different assumptions about [a/Fe] and the Reimers 
(1975) mass loss parameter rj, noting that we can only alter 
these parameters with the BaSTI models. 

In Figure [2] we show xl a^t the best fit age and metal- 
licity against absolute magnitude for our sample. There 
are some trends in some model sets, for instance the a- 
enhanced models perform best for more luminous (and mas- 
sive) galaxies, whereas BC03, PEGASE, Starburst99 and 
BaSTI scaled-solar models perform better at lower luminos- 
ity. 
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NGC4881 (BC05) NGC4881 (BC03 fit) 




2 4 6 8 10 12 14 10000 _ 

Age of SSP (Gyr) Observed Wavelength (A) 



Figure 1. Left panel: xi against age for BC03 models of different metallicities as fit to tlie pliotometry of NGC4881. Lower riglit panel: 
the best fit model spectrum ([Fe/H] = 0.093; Age = 8.5 Gyr) with the SDSS and 2MASS passbands overplotted (blue horizontal error 
bars). Upper right panel: shows the residuals from the fit (data - model) in units of 0.1 magnitude 
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Figure 2. Best fit xi against absolute magnitude for each model set. 



5.1 Age-metallicity degeneracy 



BaSTI and PEGASE models have sufficient metallicity res- 
olution that for most galaxies, there are two, three or even 
four metallicities for which an age solution can be found 
with xi within 1.0 of the lowest. Table 7 presents these 
solutions. They have very different ages, and show the well 
known age-metallicity degeneracy in the sense that the high- 
est metallicities give the youngest ages. 

In Figure [3] we illustrate the close correspondence of 
the models of different age and metallicity for IC3501. For 
each model set, three models plotted in different colours, 
and shifted vertically to separate the traces. The values of 
[Fe/H] = -0.204 for PEGASE and [Fe/H] = -0.25 for BaSTI 
provide the best fits, xt in each case is within 1 of the best 
fit, so the outer models, covering a metallicity range of ~ 
0.4 in [Fe/H], and an age range of a factor two, cannot be 
ruled out. 

Independent of model set, we find the age and metallic- 
ity changes to be linked by the following relation: 



Slog Age 
S[Fe/H] 



1.2 ±0.25 



Without other information to constrain metallicity or 
age, it appears that using broad band photometry from 
u to K, we can measure these to A[Fe/H] ~ 0.18 and 
Alog Age ~ 0.25. Formal errors on the fitted age may be 
less than this for those model sets with a coarser metallic- 
ity grid, but these errors mean very little. Metallicities and 
ages can be constrained better by using line indices, with 
e.g. Trager et al. (2008) quoting errors of 0.02 - 0.15 dex 
on metallicity values and 0.04 - 0.20 dex on log(age) values 
derived using line indices alone. 

The substantial errors on ages are disappointing, given 
the earlier claims that pairs of, for example, optical - near- 
infrared colours can break the age-metallicity degeneracy. 
Analysis of the colour-colour grids reveals why substantial 
errors should be expected, particularly for old stellar popula- 
tions. In the BaSTI simple stellar population grids oi B — K 
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Table 7. Age-metallicity degeneracy in PEGASE and BaSTI models 



PEGASE BaSTI (Scaled Solar) BaSTl (Alpha Enhanced) 

(Model 1) (Model 2) (Model 3) (Model 1) (Model 2) (Model 3) (Model 1) (Model 2) (Model 3) 





xl 


Age 




xl 


Age 




xl 


Age 




xl 


Age 




xl 


Age 




xl 


Age 




xl 


Age 




xl 


Age 


1*1 


xl 


Age 




1C3051 


2.20 


3.00 


-0.20 


2.76 


2.00 


0.00 


2.32 


l.fiO 


0.18 


1.97 


6.50 


-0.35 


1.47 


4.50 


-0.25 


1.92 


2.00 


0.06 


1.13 


15.00 


-0.70 


0.95 


14.00 


-0.60 


1.20 


2.75 


-0.29 


NGC431S 


4.00 


9.00 


-0.40 


2.33 


5.00 


-0.20 


3.87 


2.00 


0.18 


1.59 


8.00 


-0.25 


1.76 


3.25 


0.06 


2.01 


2.25 


0.26 


2.29 


13.00 


-0.60 


2.27 


5.50 


-0.29 


3.23 


2.75 


-0.09 


NGC4515 


3.16 


7.00 


-0.40 


1.93 


4.00 


-0.20 


3.21 


2.50 


0.00 


2.63 


8.00 


-0.35 


2. 56 


6.00 


-0.25 


4.11 


2.25 


0.06 


3.54 


15.00 


-0.70 


2.88 


10.50 


-0.60 


4.31 


3.25 


-0.29 


NGC4551 


2. 85 


5.00 


0.00 


2.67 


3.00 


0.18 


4.64 


2.00 


0.40 


3.15 


8.00 


0.06 


2.63 


4.00 


0.26 


2.49 


2.75 


0.40 


3.00 


13.00 


-0.29 


2.69 


8.00 


-0.09 


3.01 


5.00 


0.05 


NGC4564 


3.41 


6.00 


0.00 


2.88 


3.50 


0.18 


3.26 


2.50 


0.40 


4.50 


9.50 


0.06 


3.31 


5.00 


0.26 


2.71 


3.00 


0.40 


2. 88 


15.00 


-0.2H 


1.40 


10.00 


-0.09 


1.44 


6.50 


0.05 


NGC4867 


2.77 


6.00 


0.00 


2.67 


3.50 


0.18 


2.71 


2.50 


0.40 


4.55 


9.00 


0.06 


3.57 


5.00 


0.26 


2.82 


3.00 


0.40 


1.52 


15.00 


-0.29 


0.74 


9.50 


-0.09 


0.57 


6.00 


0.05 


NGC4872 


2.25 


6.00 


0.00 


3.03 


3.00 


0.18 


3.13 


2.00 


0.40 


3.35 


8.50 


0.06 


2.83 


4.50 


0.26 


2.38 


3.00 


0.40 


1.10 


14.50 


-0.29 


0.74 


9.00 


-0.09 


0.86 


5.50 


0.05 


NGC4S71 


2.41 


8.00 


0.00 


2.51 


3.50 


0.18 


3.38 


2.50 


0.40 


2.74 


10.00 


0.06 


2.20 


6.00 


0.26 


2.07 


3.00 


0.40 


1.63 


15.00 


-0.29 


0.86 


10.00 


-0.09 


1.24 


6.50 


0.05 


NGC4873 


2.62 


10.00 


-0.2Q 


1.86 


4.50 


0.00 


2.74 


3.00 


0.18 


2.60 


7.00 


0-06 


2.05 


3.00 


0.26 


2.90 


2.75 


0.40 


0.99 


11.00 


-0.29 


1.30 


7.00 


-0.09 


1.80 


3.75 


0.05 


NGC4473 


3.43 


8.00 


O.DO 


3.39 


3.50 


0.18 


3.81 


2.50 


0.40 


3.96 


9.50 


0.06 


2.81 


5.00 


0.26 


2.33 


3.00 


0.40 


3.74 


15.00 


-0.29 


2.39 


10.00 


-0.09 


2.59 


6.50 


0.05 


NGC48S1 


2.27 


6.00 


0.00 


2.22 


3.50 


0.18 


2.70 


2.50 


0.40 


3.77 


9.50 


0.06 


2.93 


5.00 


0.26 


2.54 


3.00 


0.40 


1.50 


15.00 


-0.29 


0.79 


10.00 


-0.09 


0.87 


6.50 


0.05 


NGC4839 


3.39 


9.00 


0.00 


2.25 


4.00 


0.18 


3.97 


3.00 


0.40 


3.95 


7.50 


0.26 


2.70 


5.00 


0.40 








1.33 


12.50 


-0.09 


1.18 


8.50 


0.05 








NGC4874 


3.25 


9.00 


0.00 


2.54 


4.50 


0.18 








3.30 


7.00 


0.26 


2.20 


5.00 


0.40 








1.31 


12.00 


-0.09 


1.47 


8.00 


0.05 








NGC4889 


4.82 


10.00 


0.00 


■A. 24: 


6.00 


0.18 


3.67 


3.00 


0.40 


G.20 


15.00 


0.06 


5.05 


9.00 


0.26 


4.37 


6.50 


0.40 


1.73 


15.00 


-0.09 


1.35 


10.50 


0.05 









Notes to Table 7: We tabulate up to 3 almost equally well fitting PEGASE, BaSTI ([a/Fe]— 0, r) — 0.4) and BaSTI ([a/Fe]— .4, 77 — 0.4) models, illustrating how models with different age and 
metallicity provide very similax quality of fits to the data. 
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Figure 3. Illustration of the age-metallicity degeneracy, we show the best age fits at four metallicities for IC3501, using PEGASE models 
(left panel) and BaSTI scaled-solar models (right panel). In each case the three older and more metal poor models are offset vertically, 
for clarity. 



vs. J — K used by James et al. (2006), the B — K colour, 
which is the principal age indicator, changes more between 
3 and 5 Gyr than it does between 5 and 14 Gyr, and over 
the range 10-14 Gyr, the change in B — K is typically only 
0.02 - 0.03 mag. In addition, for high metallicity and old 
populations, increasing age leads to redder J — K colours, 
so while the grid is not degenerate, the age and metallicity 
vectors are far from orthogonal. The same trends are found 
for theV - I YS.V - K grids of Puzia et al. (2002), which 
show a larger change in colours from 1-2 Gyr than from 5 

- 15 Gyr. Thus, while such techniques can be very sensitive 
to the presence of even a small mass fraction of young or 
intermediate-age stars, the age discrimination for old popu- 
lations is poor given typical photometric errors. Passbands 
blueward of our range may help (Kaviraj et al. 2006), as 
the models presented in Table 7 deviate from each other 
blueward of u, but uncertainties due to modelling of the 
Horizontal Branch become more serious. 

This analysis emphasises, on the other hand, the im- 
portance of having models with sufficient age and metallicity 
resolution to identify the uncertainties in the derived param- 
eters, and to investigate whether, for instance, the quality of 
the fits distinguishes between the fundamental parameters 
of the model sets, or whether the data demand more than 
an SSP to fit. 

5.2 Abundance ratios and the morphology of the 
Horizontal Branch 

The BaSTI [a/Fe] = -)-0.4 models are the only a-enhanced 
models that we investigate. For the Coma galaxies they pro- 
vide better fits than any scaled solar models, occasionally 
excepting BC03, but for the lower-luminosity Virgo galaxies 
they are less successful. The pattern of residuals shows that 
the overprediction in r is less than for the scaled-solar mod- 
els, although they still predict (r — i) to be too blue, by 0.02 

- 0.05 magnitude. However the main discrepancies are in 
{u — g), which is predicted too blue by 0.05 - 0.1 magnitude. 
The age-metallicity degeneracy still exists, with the three 
most metal rich models ([Fe/H] = -0.29, -0.09 and -1-0.05) 
showing minima with very similar values of xty but at ages 
which cover a range of a factor three, as can be seen in the 



fits to NGC4881 in the right panel of Figure |4l An extreme 
example can be seen in the left panel of this Figure, which 
shows the fits for IC3501. The metal-poor ([Fe/H]= -0.60 
and -0.70), [a/Fe] = 0.4; rj = 0.4 models fit very well at old 
ages due to a sharp transition from red to blue HB morphol- 
ogy at 8 - 14 Gyr. For other models this transition occurs 
at ages greater than the expected age of the universe. Using 
our broad-band photometry alone these models are almost 
indistinguishable from the ~ 2 Gyr, [Fe/H] ~ 0.29 models 
in this set (as Figure |4]shows), and from moderately metal- 
poor scaled-solar models. Similar effects are seen in the fits 
to NGC4318 and NGC4515. 

In general, the best fit ages for the a-enhanced models 
are greater than scaled-solar, with A{logAge) — 0.34 ±0.12, 
a mean factor of two. 

The Reimers (1975) mass-loss parameter rj has less ef- 
fect, apart from in the metal-poor a-enhanced models for the 
lowest mass galaxies discussed above. In other cases, when 
different metallicities have very close minimum xt , changing 
rj will occasionally favour a different solution which, because 
of the age-metallicity degeneracy, will then have different 
age. For the [a/Fe] = -1-0.4 models the i] — 0.2 solution can 
be one age step (500 Myr) younger at the same metallicity. 

Maraston (2005) BHB models have very high values of 
rj, between 0.45 and 1.0, and are designed to have a blue 
HB morphology even for metal-rich (0.5 - 2.0 solar) metal- 
licities. These are available for very old populations (10 - 
15 Gyr) . The ages of the best fits are constrained to be old, 
and the best fit metallicities are in the range [Fe/H] = - 
0.33 to -1-0.35 corresponding well with values derived from 
line indices. However the xt values are still high compared 
with other model sets. These models all underpredict the i 
band flux by around 0.1 magnitude, and have mixed success 
at predicting {u - g), which in the more massive galaxies is 
predicted too red by 0.1 mag. 



6 DISCUSSION 

6.1 Effect of extinction estimates 

We have applied the extinction corrections of Schlegel et 
al. (1998) but it is important to estimate the uncertainty in 
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Figure 4. Age-metallicity degeneracy in tiie a-enhanced models, 
we show the best age fits at five metaUicities to IC3501. This is 
an extreme case, with an allowed age range of a factor 10, due 
to the transition in HB morphology at ages > 8 Gyr in the high 
rj, high [o/Fe], low metallicity models. Again the older and more 
metal poor models are shifted vertically for clarity. 



the derived population parameters which results from errors 
in the extinctions. To estimate the effect of uncertainties 
in the extinction estimates, for a limited set of the stellar 
models, we have calculated the best fit age, metallicity and 
xt, assuming that Ay is 0.1 magnitude above and below 
the Schlegel et al. value. In some cases this would imply a 
negative Ay, however we use these values only to estimate 
the effect of errors in Av upon the derived parameters. 

In Table 13] we present the results for the BC03 models. 
BC03 has a coarse metallicity grid, and only in one case does 
a change of extinction correction affect the best fit metallic- 
ity. For the remaining cases, there is a clear and unsurprising 
degeneracy between the extinction correction and age. As- 
suming an additional 0.1 magnitudes above the Schlegel et 
al. value leads to a derived age a factor ~1.3 lower than the 
nominal value. However the Galactic foreground extinctions 
are not thought to be this uncertain, and Burstein & Heiles 
(1982) extinctions for our galaxies differ by less than 0.05 
mag in Av from the Schlegel et al. values, except for IC3051 
and NGC4318 where they are ~ 0.1 mag lower. Uncorrected 
internal extinction, which would strictly need to be mod- 
elled as embedded rather than as a screen, might lead to 
ages being overestimated. 



6.2 Metallicity determination 

In Figure [5] we plot the lowest value of xt, for any age, 
against [Z/H] for each model set. For BaSTI a-enhanced 
models [Z/H] = [Fe/H] + 0.35, for all other models we as- 
sume [Z/H] ~ [Fe/H]. 

With the exception of the SPEED and Maraston models 
there is broad agreement between the models sets on [Z/H]. 
However the broad minima in the Xv distributions illustrate 
further the significant uncertainties in estimates of [Z/H] or 
[Fe/H] from broad-band photometry. 



6.3 Comparison with spectroscopic ages and 
metaUicities 

A useful comparison sample for the present work is pro- 
vided by Trager et al. (2008), who present age, metallicity 
and [a/Fe] estimates from spectroscopic line index measure- 
ments for a sample of 12 Coma cluster early-type galaxies, 
including 5 from our sample. However, it is important to 
note that the models used by Trager et al. (2008) are based 
on isochrones from Worthey (1994), unlike all of the mod- 
els tested in this paper, so differences found may in part 
be due to the underlying stellar models. The quantitative 
effects on model predictions of line indices and spectral en- 
ergy distributions due to changes in models are uncertain. 
Effects resulting from e.g. differing temperature and metal- 
licity scales are being studied by Percival & Salaris (2009, 
in preparation). 

In table [5] we reproduce the age and metallicity mea- 
surements from Table 5 of Trager et al. (2008) . The quanti- 
ties they tabulate are [Z/H] and [E/Fe], where E stands for 
Enhanced elements. [E/Fe] does not precisely correspond to 
[a/Fe] in the BaSTI models, as the elements enhanced are 
somewhat different, as discussed by Trager et al. (2000a). 
However the key element for measuring [E/Fe] is Mg, which 
is an a-element. 

Trager et al. [Z/H] values are consistent with [Fe/H] 
as found from BaSTI and PEGASE scaled solar models, 
and, using the relation [Z/H] = [Fe/H]-|-0.35 for the BaSTI a- 
enhanced models, they are consistent with these also. BC03 
do not have sufficient metallicity resolution for an adequate 
comparison. Trager et al. ages are older than the best fits for 
BaSTI and PEGASE scaled solar models, but younger than 
those found for the a-enhanced models. This might imply 
that models with [a/Fe] ~ -1-0.2 could provide both a better 
fit to the data and a better match to Trager et al. (2008), 
but this is impossible to test until such models are available. 

Trager et al. (2008), and indeed much of the current 
work on stellar populations using spectroscopic data, uses 
the Lick index system. This is defined from fairly low resolu- 
tion spectra (~8.4 A FWHM), however for massive galaxies 
the effective resolution of the spectra is limited to this by the 
internal velocity dispersion (cr ~ 200km/s). Ages measured 
on the Lick system have been shown to suffer from age- 
metallicity degeneracy, and also to be somewhat dependent 
upon which Balmer line index is used (e.g. Puzia et al. 2005; 
Brodie et al. 2005). Alternative index definitions have been 
proposed (Vazdekis & Arimoto 1999; Cervantes & Vazdekis 
2009) which to some extent lift this degeneracy. In a fu- 
ture paper, we will compare the full, high-resolution BaSTI 
model spectra (Percival et al. 2009), with higher resolution 
spectra of lower velocity dispersion galaxies (e.g. Smith et 
al. 2009) to derive single or composite stellar population pa- 
rameters. 



6.4 Comparison with predictions from SBF 
magnitudes 

Lee et al. (2009), using BaSTI isochrones, investigate predic- 
tions for the 7-band Surface Brightness Fluctuation (SBF) 
magnitudes from scaled solar and a-enhanced isochrones. 
Using data from Blakeslee et al. (2009), they find the a- 
enhanced isochrones to be the better fit to the redder and 
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Table 4. Effect of different extinction estimates on the derived parameters for BC03 models 



Ay Schlegel-0.1 Schlegel extinction Ay Schlegel+0.1 

Xl Age [^] xl Age xl Age [^] 
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NGC4515 


1.80 


3.00 


0.093 


2.06 


2.75 


0.093 


2.86 


2.25 


0.093 


NGC4551 


1.74 


9.00 


0.093 


1.56 


7.50 


0.093 


1.39 


5.50 


0.093 


NGC4564 


2.54 


11.00 


0.093 


1.89 


8.50 


0.093 


1.78 


7.00 


0.093 


NGC4867 


2.30 


2.75 


0.560 


2.12 


8.00 


0.093 


1.88 


6.00 


0.093 


NGC4872 


1.62 


10.00 


0.093 


1.23 


8.00 


0.093 


0.93 


6.00 


0.093 


NGC4871 


1.20 


11.00 


0.093 


0.94 


8.50 


0.093 


1.01 


7.00 


0.093 


NGC4873 


1.19 


8.50 


0.093 


0.83 


6.50 


0.093 


0.43 


4.50 


0.093 


NGC4473 


1.87 


11.00 


0.093 


1.53 


8.50 


0.093 


1.59 


7.00 


0.093 


NGC4881 


1.56 


10.50 


0.093 


1.13 


8.50 


0.093 


1.05 


7.00 


0.093 


NGC4839 


2.55 


14.50 


0.093 


1.76 


10.00 


0.093 


1.56 


8.00 


0.093 


NGC4874 


1.89 


13.50 


0.093 


1.20 


9.50 


0.093 


1.20 


8.00 


0.093 


NGC4889 


3.55 


3.25 


0.560 


2.25 


12.00 


0.093 


1.78 


9.00 


0.093 



Notes to Table ^ Effect of extinction estimate upon the results for BC03 models. We present the derived parameters under the 
assumptions of the Schlegel et al. extinctions presented in Table 1 and adopted throughout this paper, and also for Ay values 0.10 

magnitudes above and below these. Ages are in Gyr. 



Table 5. Comparison with spectroscopic measurements from Trager et al. (2008) 





Age (Gyr) 




Trager 
[Z/H] 


et al. 

<^\Z/H] 


[E/Fe] 


0"[E/Fe] 


NGC4867 


3.0 


0.3 


0.54 


0.04 


0.20 


0.01 


NGC4872 


4.8 


0.4 


0.36 


0.02 


0.18 


0.01 


NGC4871 


4.5 


0.4 


0.36 


0.04 


0.14 


0.01 


NGC4873 


4.5 


0.4 


0.32 


0.04 


0.19 


0.01 


NGC4874 


7.9 


1.0 


0.38 


0.04 


0.17 


0.01 



Notes to Table\5jj Where Trager et al. (2008) give asymmetric error bars, we quote here the larger value. 



more massive galaxies, whereas the scaled-solar models work 
better for bluer and less massive galaxies. They attribute 
this to the effect of oxygen enhancement on the upper RGB 
and AGB. In general this result agrees with what we see 
in our data, for instance in Table [3] we see that the a- 
enhanced models provide lower values of for the most 
massive galaxies (this table is ordered by increasing luminos- 
ity). However there is still some work required to reconcile 
models, SBF magnitudes, and integrated colours, and indeed 
Lee et al. show that there are significant differences between 
Padova and BaSTI model predictions of (V^ — /) and /-band 
SBF magnitude, for the same scaled-solar abundances. 



6.5 Effect of the IMF 

In this study we have examined the effect of the choice of 
model set on the derived SSP parameters, however we have 
not examined different IMFs, having used the Kroupa (2001) 
IMF for all but one of our datasets, nor have we examined in 
detail the parametrisation of the treatment of the TP- AGB. 
Conroy et al. (2009) examined the effect of the TP- AGB 
and IMF using a dataset containing a much larger sample, 
but at lower photometric precision. It would be instructive 
to analyse our data with the Conroy et al. models, further- 
more in a future paper we shall investigate, using the BaSTI 
models, the effect of different parametrisations of the IMF 
upon the derived colours. 



7 CONCLUSIONS 

We find that SSP models provide good fits to broad-band 
photometry from u to Ks bands, for a variety of morphologi- 
cally early-type galaxies. Although we have tested only SSP 
models, we can at least say that the broad-band data do 
not require a more complicated star formation history. For 
galaxies brighter than My ~ —20.5 the best fits are usually 
provided by fully self-consistent a-enhanced models, which 
are currently only available for BaSTI. This is consistent 
with a scenario in which the stars in massive galaxies formed 
over a relatively short time interval (Thomas et al. 2005). 
For galaxies fainter than My ~ —20.5 although the popu- 
lar Bruzual & Chariot (2003) models provide in most cases 
objectively the best fit to the data, other scaled-solar abun- 
dance models including PEGASE, GALEV, Starburst99 and 
BaSTI provide very similar solutions of very similar good- 
ness of fit. This is consistent with formation of the stellar 
content of these galaxies over a longer time interval. We 
find that SPEED and Maraston (2005) models do not pro- 
vide good fits to the data, and that the best fits they do 
provide occur at unrealistically low metal abundance and 
unrealistically young age respectively. 

Broad-band photometry from it to Ka, in the presence 
of realistic photometric errors, does not fully break the age- 
metallicity degeneracy, and the models which have sufficient 
metallicity resolution show that there are strongly anticor- 
related uncertainties of ~ 0.25 in log age, and ~ 0.18 in log 



© 2008 RAS, MNRAS OOO.nifTTl 



Colours and stellar populations of galaxies 15 



IC3501 



10 



+ + BoSlf ' 

X * BC03 

o o M05_rhb 

A A M05_bhb 

□ □ BoSTLoe 

X X Pegose 

» i S99 

+ + SPEED 

>K )KGALEV 



■■. .61 



■1.0 



10 



-0.5 
[Z/H] 
NGC4881 



0.0 



0.5 



+ + BoSTI 

X X BC03 

o o M05_rhb 

A A M05_bhb 

□ □ BaSTLoe 

X X Pegose 

A 4 S99 

+ -I-SPEED 

)K )KGALEV 

J , 1 I , 1 1 I. 



-1.0 



-0.5 



0.0 



0.5 



[Z/H] 



NGC4515 



10 r-f-. 



+ + BoSTI 

X X BC03 

9 o M05_rhb 

A A M05_bhb 

□ □ BoSTLoe 

X X Pegose 

» i S99 

+ + SPEED 

>K )KGALEV 



-1 .0 



-0.5 
[Z/H] 

NGC4874 



0.0 



0.5 



10 



-+'\ 



+ + BoSTI 

X ^ BC03 

o o M05_rhb 

A A M05_bhb 

□ □ BoSTLoe 

X X Pegose 

i 1 S99 

+ -I-SPEED 

>K )KGALEV 

J 1 1 I , 1 1 ^ 



-1 .0 



-0.5 
[Z/H] 



0.0 0.5 



Figure 5. x5 for best fitting age against [Z/H] for eacli model set, for IC3051 (top left), NGC4515 (top right), NGC4881 (bottom left) 
and NGC4874 (bottom right). Models are denoted by the symbols, and diff'erent metallicities connected by the line colours and styles, 
as denoted in the key in the box. 



metallicity, on the derived population parameters. BaSTI 
Q-enhanced models fit at older ages than the scaled-solar 
models, but age-metallicity degeneracy is still present. 

BaSTI models use stellar tracks with two values of both 
[a/Fe] and the Reimers (1975) mass-loss parameter rj, which 
controls horizontal branch morphology. We find that the lat- 
ter has a much smaller effect than the former, but the range 
of this parameter that we have explored in this paper is 
small. 

To improve our knowledge of stellar population pa- 
rameters, it is important to investigate complementary 
techniques. These will include extending the wavelength 
range used further into the ultraviolet with HST/WFC3, 
which will provide an additional tool for breaking the age- 
metallicity degeneracy. However to utilize the UV photom- 
etry it is particularly important to understand the contri- 
bution of the regions of the HR diagram which contribute 
most to the UV flux, in particular the Horizontal Branch, 
Blue Stragglers, and any contribution from a young stellar 
component, for instance from merger induced star forma- 



tion. Further complementary techniques include studies of 
the properties of star clusters associated with the galaxies, 
which sometimes show subpopulations indicating different 
stellar population parameters, and study of surface bright- 
ness fluctuations at near infra-red wavelengths, which are a 
powerful probe of the structure of the AGB and upper RGB. 
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